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We report results from a blind analysis of the final data taken with the Cryogenic Dark Matter 
Search experiment (CDMS 11) at the Soudan Underground Laboratory, Minnesota, USA. A total 
raw exposure of 612 kg-days was analyzed for this work. We observed two events in the signal region; 
based on our background estimate, the probability of observing two or more background events is 
23%. These data set an upper limit on the Weakly Interacting Massive Particle (WIMP)-nucleon 
elastic-scattering spin-independent cross-section of 7.0 x 10"'*'' cm'^ for a WIMP of mass 70GeV/c'^ 
at the 90% confidence level. Combining this result with all previous CDMS II data gives an upper 
limit on the WIMP-nucleon spin-independent cross-section of 3.8 x 10"^'* cm^ for a WIMP of mass 
70 GeV/c^ We also exclude new parameter space in recently proposed inelastic dark matter models. 

PACS numbers: 14.80.Ly, 95.35. -|-d, 95.30.Cq, 95.30.-k, 85.25.0j, 29.40.Wk 



Cosmological observations [T] have led to a concor- 
dance model of the universe where ~85% of matter is 
non-baryonic, non- luminous and non-relativistic at the 
time of structure formation. Weakly Interacting Mas- 
sive Particles (WIMPs) ^ are a class of candidates for 
this dark matter which are particularly well motivated 
by proposed extensions to the Standard Model of parti- 
cle physics and by thermal production models for dark 
matter in the early universe [31 H |5J [5]. WIMPs, dis- 
tributed in a halo surrounding our galaxy, would coher- 
ently scatter off nuclei in terrestrial detectors [Tj [HI [5] 
with a mean recoil energy of tens of keV, presently 
limited by observation to a rate of less than 0.1 event 



kg~*day~* [5J[H1[TU]. Direct search experiments seek re- 
coil signatures of these interactions and have achieved the 
sensitivity to begin testing the most interesting classes of 
WIMP models [HI HH [H [U . 

The Cryogenic Dark Matter Search (CDMS II) exper- 
iment, located at the Soudan Underground Laboratory, 
uses 19 Ge (^^230 g) and 11 Si (~100 g) particle detectors 
operated at cryogenic temperatures (< 50 mK) [TTl [T5] . 
Each detector is a disk ~10 mm thick and 76 mm in diam- 
eter. Particle interactions in the detectors deposit energy 
in the form of phonons and ionization. Phonon sensors 
on the top of each detector are connected to four phonon 
readout channels to allow measurement of the recoil en- 



2 



ergy and position of an event. These phonon sensors 
are also the ground reference for the ionization measure- 
ment. The electric field for the ionization measurement 
is formed by applying a voltage bias to the bottom detec- 
tor surface, which is segmented into two concentric elec- 
trodes. Events having an ionization signal in the outer 
ionization channel of the detector are excluded, defining 
an ionization fiducial volume. The detectors are grouped 
into five towers, each tower containing six detectors. De- 
tectors are identified by their tower number (T1-T5) and 
their position within that tower (Z1-Z6). A direct line of 
sight between adjacent detectors in a tower allows iden- 
tification of events scattering between detectors. 

The ratio of the ionization to recoil energy ("ioniza- 
tion yield") provides event-by-event rejection of electron- 
recoil events to better than lO^"' misidentification. Es- 
sentially all of the misidentified electron recoils are "sur- 
face events" occurring within the first few /im of the de- 
tector surface, which can suffer from sufficiently reduced 
ionization collection to be misclassified as nuclear recoils. 
Due to interactions of phonons in the surface metal lay- 
ers, these surface events have faster-rising phonon pulses 
than events occurring within the bulk of the detectors 
("bulk events"). Hence, we use phonon pulse timing pa- 
rameters to improve rejection of surface events to ob- 
tain an overall (yield and timing) misidentification prob- 
ability of less than 10^^ for electron recoils. To atten- 
uate external environmental radioactivity and to reject 
events caused by cosmogenic muons, the detectors are 
surrounded by layers of lead and polyethylene shielding 
and an active muon veto |15j . 

Data taken during four periods of stable operation be- 
tween July 2007 and September 2008 were analyzed for 
this work. Between each data period the cryostat was 
warmed up for maintenance of the cryogenic system. All 
30 detectors were used to identify particle interactions, 
but only Ge detectors were used to search for WIMP 
scatters. Five Ge detectors were not used for WIMP de- 
tection because of poor performance or insufficient cali- 
bration data; four more detectors were similarly excluded 
during subsets of the four periods. We excluded Si de- 
tectors in this analysis due to their lower sensitivity to 
coherent nuclear elastic scattering. 

A subset of events were analyzed to monitor de- 
tector stability and identify periods of poor detector 
performance. Data quality criteria were developed on 
a detector-by-detector basis using Kolmogorov-Smirnov 
tests performed on parameter distributions. Our detec- 
tors require regular neutralization |15j to maintain full 
ionization collection. We monitor the yield distribution 
and remove periods with poor ionization collection. Af- 
ter these data quality selections, the total exposure to 
WIMPs considered for this work was 612 kg-days. 

The data selection criteria ("cuts") that define the 
WIMP acceptance region were calculated from calibra- 
tion sets of electron and nuclear recoils obtained during 



regular exposures of the detectors to 7-ray (^^'^Ba) and 
neutron (^^^Cf ) sources. A 356 keV 7-line from the ^^^^Ba 
source was used to calibrate the ionization and phonon 
energy scales of each detector. Nuclear recoils from the 
elastic scattering of neutrons emitted by the ^^^Cf source 
were used to define the WIMP signal region, taken to be 
the 2(7 band about the mean neutron ionization yield in 
the yield vs recoil energy plane (Fig. [T]) . Neutron activa- 
tion of ^°Ge results in the emission of 10.36 keV x-rays. 
The width of this x-ray line determined the ionization 
energy resolution at the lOkeV analysis threshold to be 
< 0.4 keV on all detectors considered. 

Phonon pulse shapes vary with event position and en- 
ergy in our detectors. To enhance surface-event rejection, 
we characterized the variation of the resulting timing and 
energy estimators using our electron-recoil calibrations 
and corrected for it in both WIMP-search and calibra- 
tion data. To ensure this correction works well at large 
radius, where position reconstruction degeneracies can 
lead to miscorrection, we included calibration electron- 
recoil events outside the ionization fiducial volume in the 
correction calculation. This technique obviated the addi- 
tional fiducial volume cut in phonon position reconstruc- 
tion parameters used in past analyses [11] . 

Figure [1] demonstrates our surface-event rejection cut. 
We optimized this cut using calibration sets of nuclear 
and surface electron recoils from the ^^^Cf and ^'^^Ba ex- 
posures. We verified that the sum of the rise time of the 
largest phonon pulse plus its delay relative to the ion- 
ization signal provided the best discrimination between 
surface events and nuclear recoils. Surface-event rejec- 
tion criteria based on this discriminator were tuned on 
the calibration data by maximizing the expected sensi- 
tivity for a 60GeV/c2 WIMP. 

A blind analysis was performed, in which cuts were de- 
veloped without looking at events that might appear in 
or near the signal region. Candidate WIMP-scatters were 
required to be within 2a of the mean ionization yield of 
nuclear recoils and at least 3(t away from the mean ioniza- 
tion yield of electron recoils, have recoil energy between 
10 and 100 keV, and have ionization energy at least 4.5ct 
above the noise. They are required to occur within the 
detector fiducial volume, satisfy data quality criteria and 
pass the surface-event rejection cut. Since WIMPs are 
expected to interact only once ( "single-scatter event" ) in 
the experimental apparatus, a candidate event was fur- 
ther required to have energy deposition consistent with 
noise in 29 detectors. Additionally, we required the ab- 
sence of significant activity in the surrounding scintillator 
veto shield during a 200-/is window around the trigger. 

The efficiency of our analysis cuts for nuclear recoils 
was measured as a function of energy using both neutron- 
calibration and WIMP-search data. The fiducial vol- 
ume estimate is corrected for neutron multiple-scattering 
based on simulations. Our efficiency for signal events has 
a maximum of -32% at 20keV. It falls to -25% both 
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FIG. 1: The power of the primary background discrimination 
parameters, ionization yield and phonon timing, is illustrated 
for a typical detector using in situ calibration sources. Shown 
are bulk electron recoils (red points), surface electron events 
(black crosses) and nuclear recoils (blue circles) with recoil 
energy between 10 and 100 keV. Top: Ionization yield ver- 
sus recoil energy. The solid black lines define bands that are 
2a from the mean electron- and nuclear-recoil yields. The 
sloping magenta line indicates the ionization energy thresh- 
old while the vertical dashed line is the recoil energy analy- 
sis threshold. The region enclosed by the black dotted lines 
defines the sample of events that are used to develop surface- 
event cuts. Bottom: Normalized ionization yield (number of 
standard deviations from mean of nuclear recoil band) versus 
normalized timing parameter (timing relative to acceptance 
region) is shown for the same data. Events to the right of 
the vertical red dashed line pass the surface-event rejection 
cut for this detector. The solid red box is the WIMP signal 
region. (Color online.) 

at 10 keV, due to ionization threshold and flaring of the 
electron-recoil band; and at 100 keV, due to a drop in 
fiducial volume. The spectrum- averaged equivalent ex- 
posure for a WIMP of mass 60GeV/c^ is 194.1 kg-days. 

Neutrons with energies of several MeV can generate 
single-scatter nuclear recoils that are indistinguishable 
from possible dark matter interactions. Sources of neu- 
tron background include cosmic-ray muons interacting 
near the experimental apparatus (outside the veto), ra- 
dioactive contamination of materials, and environmen- 



tal radioactivity. We performed Monte Carlo simula- 
tions of the muon-induced particle showers and sub- 
sequent neutron production with Geant4 [lH [17] and 
FLUKA [T51 [12]. The cosmogenic background is esti- 
mated by multiplying the observed number of vetoed sin- 
gle nuclear recoils in the data by the ratio of unvetoed 
to vetoed events as determined by cosmogenic simula- 
tion. This technique predicts 0.04l;Q Q3(stat) events in 
this WIMP-search exposure. 

Samples of our shielding and detector materials were 
screened for U and Th daughters using high purity ger- 
manium 7 counters. In addition, a global 7-ray Monte 
Carlo was performed and compared to the electromag- 
netic spectrum measured by our detectors. The contam- 
ination levels thus determined were used as input to a 
Geant4 simulation to calculate the number of neutrons 
produced from spontaneous fission and (a, n) processes, 
assuming secular equilibrium. The estimated background 
is between 0.03 and 0.06 events and is dominated by U 
spontaneous fission in the copper cans of the cryostat. 
The radiogenic neutron background originating from the 
surrounding rock is estimated to be negligibly small com- 
pared to other sources. 

The number of misidentified surface events was esti- 
mated by multiplying the observed number of single- 
scatter events failing the timing cut inside the 2a nuclear- 
recoil band with the ratio of events expected to pass the 
timing cut to those failing it ( "pass- fail ratio" ) . The for- 
mer was estimated using observed counts from a previous 
analysis [11] , and the latter was estimated using three dif- 
ferent methods. The first method computed the pass- fail 
ratio from events that reside within the 2a nuclear-recoil 
band and multiply scatter in vertically adjacent detec- 
tors ( "multiple scatter events" ) . The second method esti- 
mated the pass-fail ratio from multiple-scatter events sur- 
rounding the nuclear-recoil band ("wide-band events"). 
Wide-band events have different distributions in energy 
and in detector face (ionization- or phonon- side) from 
nuclear-recoil band events, affecting the pass-fail ratio. 
To account for these differences, the pass-fail ratio of 
these events was corrected using the face and energy dis- 
tributions of events observed in the nuclear-recoil band 
that failed the timing cut. A third, independent estimate 
of the pass-fail ratio was made using low-yield, multiple- 
scatter events in ^■^•^Ba calibration data, again adjusted 
for differences in energy and detector-face distributions. 
All three estimates were consistent with each other and 
were thus combined to obtain an estimate prior to un- 
blinding of 0.6 ± O.l(stat) surface events misidentified as 
nuclear recoils. 

Upon unblinding, we observed two events in the WIMP 
acceptance region at recoil energies of 12.3 keV and 
15.5 keV. These events are shown in Figs. [2] and [3) 

The candidate events occurred during periods of nearly 
ideal experimental performance, are separated in time 
by several months, and occur in different towers. How- 
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FIG. 2: Ionization yield versus recoil energy for events pass- 
ing all cuts, excluding yield and timing. The top (bottom) 
plot shows events for detector T1Z5(T3Z4). The solid red 
lines indicate the 2cr electron and nuclear recoil bands. The 
vertical dashed line represents the recoil energy threshold and 
the sloping magenta dashed line is the ionization threshold. 
Events that pass the timing cut are shown with round mark- 
ers. The candidate events are the round markers inside the 
nuclear-recoil bands. (Color online.) 



FIG. 3: Normalized ionization yield (number of standard de- 
viations from mean of nuclear recoil band) versus normalized 
timing parameter (timing relative to acceptance region) for 
events passing all cuts, excluding yield and timing. The top 
(bottom) plot shows events for detector T1Z5(T3Z4). Events 
that pass the phonon timing cut are shown with round mark- 
ers. The solid red box indicates the signal region for that 
detector. The candidate events are the round markers inside 
the signal regions. (Color online.) 



ever, a detailed study revealed that an approximation 
made during the ionization pulse reconstruction degrades 
the timing-cut rejection of a small fraction of surface 
events with ionization energy below ^--^6 keV. The can- 
didate event in T3Z4 shows this effect. Such events 
are more prevalent in WIMP-search data than in the 
data sets used to generate the pre-blinding estimate of 
misidentified surface events. A refined calculation, which 
accounts for this reconstruction degradation, produced a 
revised surface-event estimate of 0.8±0.1(stat)±0.2(syst) 
events. The systematic uncertainty is dominated by our 
assumption that the pass-fail ratio for multiple scatter 
events is the same as that for single scatter events. Based 
on this revised estimate, the probability to have observed 
two or more surface events in this exposure is 20%; in- 
clusion of the neutron background estimate increases this 
probability to 23%. These expectations indicate that the 
results of this analysis cannot be interpreted as significant 
evidence for WIMP interactions, but we cannot reject ei- 
ther event as signal. 

To quantify the proximity of these events to the 
surface-event rejection threshold, we varied the timing 



cut threshold of the analysis. Reducing the revised ex- 
pected surface-event background to 0.4 events would re- 
move both candidates while reducing the WIMP expo- 
sure by 28%. No additional events would be added to 
the signal region until we increased the revised estimate 
of the expected surface-event background to 1.7 events. 

We calculate an upper limit on the WIMP-nucleon 
elastic scattering cross-section based on standard galactic 
halo assumptions |10j and in the presence of two events 
at the observed energies. We use the Optimum Interval 
Method [22] with no background subtraction. The result- 
ing limit shown in Fig. |4]has a minimum cross section of 
7.0 X lO"^"* cm^ (3.8 x 10"^'' cm^ when combined with 
our previous results) for a WIMP of mass 70GeV/c^ . 
The abrupt feature near the minimum of the new limit 
curve is a consequence of a threshold-crossing at which 
intervals containing one event enter into the optimum in- 
terval computation |22j . An improved estimate of our 
detector masses was used for the exposure calculation of 
the present work; a similar correction (resulting in a ~9% 
decrease in exposure) was applied to our previous CDMS 
result [TT] shown in Fig. |4] While this work represents 
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FIG. 4: 90% C.L. upper limits on the WIMP-nucleon spin- 
independent cross section as a function of WIMP mass. The 
red (upper) solid line shows the limit obtained from the ex- 
posure analyzed in this work. The solid black line shows 
the combined limit for the full data set recorded at Soudan. 
The dotted line indicates the expected sensitivity for this ex- 
posure based on our estimated background combined with 
the observed sensitivity of past Soudan data. Prior results 
from CDMS [11], EDELWEISS II [S], XENONIO [13], and 
ZEPLIN III are shown for comparison. The shaded re- 
gions indicate allowed parameter space calculated from cer- 
tain Minimal Supersymmetric Models [201 121| (Color online.) 

a doubling of previously analyzed exposure, the observa- 
tion of two events leaves the combined limit, shown in 
Fig. [4] nearly unchanged below 60GeV/c^ and allows 
for a modest strengthening in the limit above this mass. 

We have also analyzed our data under the hypothesis 
of WIMP inelastic scattering 23J, which has been pro- 
posed to explain the DAMA/LIBRA data . We com- 
puted DAMA/LIBRA regions allowed at the 90% C.L. 
following the ^ goodness-of-fit technique described in 
(25] . without including channeling effects [26]. Limits 
from our data and that of XENONIO [57] were com- 
puted using the Optimum Interval Method [25]. Re- 
gions excluded by CDMS and XENONIO were defined 
by demanding the 90% C.L. upper limit to completely 
rule out the DAMA/LIBRA allowed cross section in- 
tervals for allowed WIMP masses and mass splittings. 
The results are shown in Fig. [5j The CDMS data dis- 
favor all but a narrow region of the parameter space al- 
lowed by DAMA/LIBRA that resides at a WIMP mass 
of ~100GeV/c2 and mass splittings of 80-140 keV. 

The data presented in this work constitute the final 
data runs of the CDMS II experiment and double the 
analyzed exposure of CDMS II. We observed two can- 



FIG. 5: The shaded green region represents WIMP masses 
and mass splittings for which there exists a cross section com- 
patible with the DAMA/LIBRA ^ modulation spectrum 
at 90% C.L. under the inelastic dark matter interpretation 
[23 . Excluded regions for CDMS II (solid-black hatched) and 
XENONIO 27 (red-dashed hatched) were calculated in this 
work using the Optimum Interval Method. (Color online.) 

didate events. These data, combined with our previous 
results, produce the strongest limit on spin-independent 
WIMP-induced nuclear scattering for WIMP masses 
above 42GeV/c^ ruling out new parameter space. 
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